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Abstract. 

Starting from Rodrigues formula we present a general construction of raising and 
lowering operators for orthogonal polynomials of continuous and discrete variable 
on uniform lattice. In order to have these operators mutually adjoint we introduce 
orthonormal functions with respect to the scalar product of unit weight. Using the 
Infeld-HuU factorization method, we generate from the raising and lowering operators 
the second order self-adjoint differential/difference operator of hypergeometric type. 
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1. Introduction 

The factorization metliod lias become a very powerful tool to solve second order 
differential equations and its application to physical models with orthonormal basis, 
generated by creation and annihilation operators. A classical paper by Infeld and Hull 
[1] defined the method and applied it to a large class of second order Hamiltonian 
that generalizes the well known description of the non-relativistic oscillator by means 
of creation and annihilation operators. Miller [2] enlarged this method to difference 
equation and made connection to orthogonal polynomial of discrete variable. An analysis 
of the factorization types led Miller to the idea that this method is a particular case of 
the representation theory of Lie algebras. 

The two volumes of Nikiforov, and colaborators [3] [4] about classical orthogonal 
polynomials of continuous and discrete variable opened the way to a more rigurous and 
systematic approach to the factorization method. In fact Atakishiev and colaborators 
[5] [6] [7] [8] explored the application of Kravchuk, Meixner and Charlier polynomials 
to the eigenvalue problem of some dynamical system where the energy eigenvalues are 
equally spaced. This particular situation makes possible to determine the generators of 
the dynamical symmetry group. Also Smirnov [9] has used the properties of difference 
equations of hypergeometric type given in [3] to construct raising and lowering operators 
that generate orthonormal functions corresponding to Hamiltonian of different levels. 

Bangerezako and Magnus have developed the method of the factorization of 
difference operators of hypergeometric type [10] [11] [12]. They proposed two different 
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approaches for this factorization. 1) For a given operator find raising and lowering 
operators that generate a complete set of polynomials cigcnfunctions. 2) Generate from 
a factorization chain an operator having a complete set of polynomial eigenf unctions. 

We have presented two papers [13] [14] related to the construction of the creation 
and annihilation operators for the orthogonal polynomials (or functions) of continuous 
and discrete variables. The motivation for these papers was the construction of 
mathematical model for Quantum systems on discrete space-time (as the harmonic 
oscillator, the hydrogen atom, the Dirac equation) [15] and to make connection with 
standard quantum mechanics by the continuous limit. 

In this paper we follow the second approach to the factorization method of 
Bangerezako and Magnus explained before. Starting from the raising and lowering 
operators we generate the second order differential/difference equation corresponding 
to the hypergeometric functions of continuous and discrete variable. Our procedure is 
completely general and valid for all functions of this type. 

In section 2 we use the results of Nikiforov, Suslov and Uvarov [3] [4] connecting first 
order derivatives and orthogonal polynomials (as a consequence of Rodrigues' formula) 
to construct raising and lowering operators (the last one with the help of recurrence 
relations). In general these operators are not mutually adjoint with respect to the 
standard scalar product. For this reason, we introduce in section 3 the orthonormalized 
functions of hypergeometric type and then the corresponding raising and lowering 
operators are always mutually adjoint. 

In section 4 and 5 we repeat the same systematic procedure, derived from 
Rodrigues's formula, to calculate the raising and lowering operators for orthogonal 
polynomials and functions of discrete variable. It can be proved that these operators 
are also mutually adjoint. 

In section 6 we introduce the factorization method to generate the second order 
differential operator of the Sturm-Liouville type having a complete set of polynomials 
cigcnfunctions. The factorization of the raising and lowering operators fulfills (up to a 
factor) the defining equations of the Infeld-Hull method [1]. 

In section 7 we apply the same technique to the hypergeometric functions of discrete 
variable. As usual all these functions transform in the limit into the corresponding 
hypergeometric functions of continuous variable. 

It is important to make clear that the raising and lowering operators, introduced 
in sections 2 to 5, are defined with respect to one index only, namely, the degree of 
the orthogonal polynomials or the degree of the corresponding orthonormal functions. 
The same definition has been used by Atakishiyev and colaborators [6], [7], [8], by 
Bangerezako and Magnus [10], [11], [12] and by Infeld, Hull and Miller [1], [2]. Physically 
this situation corresponds in the case of quantum oscillators to the creation and 
annihilation operators with respect to the index that distingues different eingenvectors 
of the energy operators. 

With respect to the factorization techniques in the case of difference equations of 
section 7 two types of factorization can be considered [12]. Writting a hnear difference 



3 



equation of second order in the form 



d 



H{x)y{x) = M^)EM 



i=—d 



where E^. [/(^c)] = f{x + i), o? £ i E Z and Ai[x) are some scalar functions in the 
first type of factorization consists in factorizing exactly the operators H{x) + C, with 
C some constant, and the raising and lowering operators satisfying a quasi-periodicity 
condition (Spiridinov-Vinet-Zhedanov type) [19]. The second factorization technique 
consists in factorizing the operatorE^^ o [H{x) + C] with some raising and lowering 
operators that are shape- invariant (Infeld- Hull- Miller type) [1], [2]. In section 7 we 
have used the first type of factorization, but in section 6 obviously we have used the 
Infeld-HuU-Miller technique for differential equations of hypergeometric type. 

2. Raising and lowering operators for orthogonal polynomicils of continuous 
veiriable 

A polynomial of hypergeometric type yn{s) of continuous variable s satisfies two 
fundamental equations; from which one derives raising and lowering operator 

i) Differential equation 



where a{s) and r(s) arc polynomials of, at most, second and first degree respectively, 
and A„ is a constant, related to the above functions 



The differential equation can be written in the form of an eigenvalue equation of 
Sturm-Liouville type: 



where p{s) id the weight function, satisfying [a{s)p{s))' — t{s)p{s) . 

The solutions of the differential equation are polynomials that satisfy an 
orthogonality relation with respect to the scalar product 



where (i„ is a normalization constant. 

The differential equation (CI) defines an operator that is self adjoint with respect 
to this scalar product 

ii) Three term recurrence relations: 

(C2) Syn{s) = anVn+lis) + PnVnis) + InVn-lis) 

where an, (3n, In a-re constants. 



(CI) 



a{s)y'; + T{s)y'n{s) + X^ynis) = 




{a{s)p{s)y'n{s))' + Xnp{s)yn{s) = 




iii) Raising operator: From the Rodrigues formula (which is a consequence of the 
differential equation (CI) one derives a relation for the first derivative of polynomials 
yn{s) in terms of the polynomials themselves 



rn{s)yn{s) - --^yn+i{s) 



where 



T„(s) = t(s) + na'{s) 
r;(.)=r' + na": 



~2n + l 

We can modify the last equation in a more suitable form. From 
an^B^Yl(T' + ^{n + k-l)a"^ , ao = Bo 

k=0 ^ ^ 

we can prove the following identity 



an — 



an Bn t' + ^a" D„ A„ 2n 2n + 1 



ttn+i (r' + ^a") (r' + na") B^+i n X2„ Asn+i 

from which we finally get 

(C3) + — ^4^2/n(s) - a{s)y'^{s) = ^a^yn+iis) 
n T„ In 

The left side of this equation can be considered the differential operator which, 
when applied to yn{s) , gives a polynomial of higher degree. 

iv) Lowering operator: 

Introducing (C2) in (C3) we get a differential operator which, when applied to an 
orthogonal polynomial of some degree, gives another polynomial of lower degree. 

(C4) - + ^{s - Pn) + a{s)y'^{s) = ^7nyn-i{s) 

Formula (C3) and (C4) can be used to calculate solutions of the diferential equation 
(CI). In fact, if we put n = in (C4) we get yo{s) . Inserting this value in (C3) we 
obtain by iteration all the solutions of the differential operator (CI). 

The explicit expressions for orthogonal polynomials of continuous variable are given 
in Table I. The values of p{s), a{s),T{s), A„, «„, 7„, o?„ are taken from [4]. 



3. Raising and lowering operators for orthonormal functions of continuous 
variable 

From the orthogonal polynomials that satisfy a scalar product with respect to the weight 
p{s) we can construct a new functions: 

V'n(s) = d-^y^p(s}yn{s) 



5 



and obtain orthogonal functions of unit norm. Solving the last expression for ?/„(s) and 
substituing in (CI), (C2), (C3) and (C4) and using the properties of a{s) and t(s) we 
obtain the following expressions for the normalized orthogonal functions: 

i) Differential equation: 



(NCI) a(sX(s) + a'(s)V^;(s) 



1 (r(s) - a'{s) f ^ 1 _ 



4 a{s) 

which corresponds to a self-adjoint operator of Sturm- Liouville type, 
ii) Recurrence relation: 

(NC2) —anll^n+l{s) + -. ./. /-^ . 



V'n(s) + K^n{s) = 



2n dr, 



iii) Raising and lowering operators: 



(NC3) L+{s,n)^/jn{s) 

^2n _ dn+1 



n T, 



7^ + {r{s) - a'{s)) 



^n{s) - a{s)i)'^{s) 



(NC4) L-{s,n)i,n{s) 



^2n dn-1 



Putting n = in (NC4) we obtain ipo{s) , and inserting this value in (NC3) we 
obtain by iteration all the orthonormal function of hypcrgcomctric type. 

The explicit expressions for these functions are given in Table I. 

We want to make two observations. First, the operator corresponding to (NCI) 
is a self-adjoint operator of Sturm-Liouville type as can be easily checked. Secondly, 
the raising and lowering operators (NC3) and (NC4) are mutually adjoint in the case 
of Laguerre and Hermite functions. For the Jacobi and Legendre functions we have to 
multiply both operators by 2n/\2n- In fact, we have 

b 



I 



V'n+l(s) 



2n 

^2n 



L'^{s,n)ipn{s) 



ds — a, 



d 



n+l 
dn 



2n + 2 



A 



2n+2 



L {s,n + l)tjjn+i{s) 



'll)n{s)ds = -in+l-J^ 



Both integrals are equal because 7„+i = q;„-^ 
(In the case of Hermite and Laguerre functions A^/m is independent of m, for any m) 



6 



4. Raising and lowering operators for orthogonal polynomials of discrete 
variable 

A polynomial of hypergeometric type Pn{x) of discrete variable x satisfies two 
fundamental relations from which one derives raising and lowering operators. 

i) Difference equation: 

(Dl) a{x)AVPn{x) + T{x)APr,{x) + XnPnix) = 

where (t{x) and r{x) are polynomials of, at most, second and first degree, respectively. 
The forward (backward) difference operators are: 

Af{x) = fix + 1) - fix) V/(x) = fix) - fix - 1) 

This difference equation can be written in the form of an eigenvalue equation of 
Sturm Liouville type 

A [aix)pix)VPnix)] + Kpix)Pnix) = 

where pix) is a weight function satisfying 

A [a(x)p(x)] = r(,T)p(x) 
and A„ is the eigenvalue corresponding to the eigenfunction P(,t) : 

A„ = -nArix) - ^^^^^ AV(x) = -n [r' + 

The solution of the difference equation are polynomials that satisfy an orthogonality 
relation with respect to the scalar product 

fe-i 

^Pnix)Pmix)pix) = d^dnm 

x=a 

when 5mn is the Kronecker symbol and dn some normalization constant. The difference 
equation (Dl) defines an operator that is self-adjoint with respect to this scalar product. 

ii) Three term recurrence relations: 

(D2) xPnix) = anPn+iix) + (3„P„ix) + 7„P„+i(x) 
where «„, 7„ are some constants. 

iii) Raising operator 

Prom the Rodrigues formula, one derives a relation for the first difference operator 
of polynomials Pnix) in terms of the polynomials themselves. 



aix)VPr,ix) 



nrL 



rnix)Pnix) - — Pn-\-iix) 



where 



Tnix) — rix + n) + aix + n) — aix) 
Arnix) = Arix) + nA^aix) 



/ I It/ \ -^2n+l 
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because a{x) and t{x) are polynomials of at most second and first order degree 
respectively. 

We can modify the last equation to a more suitable form, as we did in the continuous 
case. From the definition 



V 1 \ 

an^ BnYlir' + -{n + k-l)a"], Uq ^ Bo 

h—n V / 



fe=0 

we have the following identity: 

2n {2n + 1) n B„ 



from which we get a more simplified version 

(D3) aix)VP^ix) = ^^Pnix) - ^Pn+iix) 

n in 

This equation defines the raising operator in terms of the backward difference. 

iv) Lowering operator 

Prom the expresion for the raising operator we can derive another lowering operator 
in terms of the forward operator. We substitute the difference operator V in (D3) for 
its equivalent V = A — VA , and then the difference equation (Dl) and the three terms 
recurrence relations (D2), with the result 



n X2n+i 2n 



Pn{x) 



(D4) {a{x) + t{x)) ^Pn{x) = 

+ ^lnPn-l{.x) 

As in the continuous cas from (D4) putting n = we get Pq{x) and inserting this 
value in (D3) we obtain by iteration all the polynomials Pn{x) satisfying (Dl). 

The explicit expressions for the orthogonal polynomials Pn{x) are given in table II. 
The values of p(a;), a{x),r{x), A„, «„, 7„, d„ are taken from [4]. 



5. Raising and lowering operators for orthonormal functions of discrete 
VEiriable 

In the last section we have a set of polynomials that arc orthogonal with respect to 
the weight function p[x) . From these polynomials we construct some functions that 
are orthogonal with respect to the unit weight, p{x) = 1 , and at the same time are 
normalized. 

Introducing this expresion in (Dl), (D2), (D3) and (D4) and using the properties 
of function a{x), t{x) and p{x), we obtain 



i) Difference equation 



(NDl) ^/{a{x)+T{x))a{x + l)(f)n{ X + 1) + a/ {a{x — 1) + t{x — 1)) (7(x)0„(x — 1) 

- {2a{x) + r{x)) ^x) + XnM^) = 
ii) Three term recursion relation: 

(ND2) ^Q;n%^0n+l(2;) + l^ln^^^(l>n-l{x) + ^(/3„ - x)(j)n{x) = 

2n dn 2m, dn 2n 



iii) Raising operator 
(ND3) L+{x,n) = 



A„rn(a;) 



n T' 



^271 dn+1 , 

In dn 



- (J[x) 



x] 



(f)n{x) + V {a{x - 1) + t{x - 1)) a{x)(f)n{x - 1) 



iv) Lowering operator 

(ND4) L-{x,n) = 



— + An + —{X - (3n) - Cr{x) - t{x) 



n tL 



4>n{x) 



+ a/ {a{x) + t{x)) a{x + l)0„(x + 1) = ^7n%^0„_i(a;) 

2n dn 

From the last two expressions we get all the solution of the difference equation 
(NDl). Putting n = in (ND4) we obtain (f)o{x) and inserting this value en (DC3) we 
obtain, by iteration, all the normalized function (j)n{x)- 

The explicit calculations for all the orthonormal functions of hypergeometric type 
are given in table II. 

As in section 3, we make two observations. Firstly, the raising and lowering 
operators (ND3) and (ND4) are mutually adjoint in the case of Krauvchuk, Meixner 
and Charlier functions. For the Hahn and Chebyshev functions we have to divide both 
by A2n/2n, therefore, they become mutually adjoint, namely. 



b-l 



'^<Pn+l{x) 



2n 



A 



2n 



L+{x,n)(f)n{x) 



b-l 

E 



2n + 2 



X 



2n+2 



L {x,n+l)(f)n+i{x) 



{x) = an 



dr,. 



Secondly, the operator corresponding to the eigenvaue A„ in (NDl) is self adjoint. 
In order to prove it, it is enough to show 

6-1 

Xl'^K^) {yR^yT7(^)M^Tl)0n(a^ + l) + V((7(x - 1) + t{x - 1)) a(x)(l)n(x - 1)} 

x=a 

b-l 

= { VH^ - 1) + - 1)) a(x)(l)i(x - 1) 

x=a 

+ y/(a(x) + t(x)) a(x + l)0i(x + 1)} 
From the orthogonality conditions a{a) = a{h) — 0, we can write 

b-l 

4>n{x) V((7(x - 1) + t{x - 1)) a{x)(t>i{x - 1) 
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6-2 



x'=a—l 
6-1 

x=a 

Similarly 

6-1 

x=a 

6-1 

= "l^nix - 1) V((7(x - 1) + r{x - 1)) a{x)(f)i{x) 



6. Factorization for differential equation of hypergeometric type 

The raising and lowering operators of section 2 and 4 will help us to factorize the 
second order differential equation of hypergeometric type into the product of two first 
order operators in agreement with the general method of Infeld and Hull [1] . 
From (NCI) we define the operator 

that satisfies H{s,n)ilJn{s) = 

We write the raising and lowering operators, (NC3) and (NC4) respectively, in the 
following way 

L+(s,n) = - cr{s)j^ 

L~{s,n) = g{s,n) + a{s)-^ 



where 



satisfying 



f{s,n-l) = g{s,n) or f{s,n)^g{s,n + l) 



which can be proved by Taylor expansion. 
Now we calculate 

L~(s, n + n) = g{s, n + n) + a{s) {/(s, n) - g{s, n + 1)}^ 
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The second term of the right side becomes zero. Substituting the values for f{s, n), 
g{s,n) and H{s,n) we get 



A 



TrAS, 



n T' 



n 



L"(s,n+ l)L+(s,n) = 

—a{s)H{s, n) 

It can be proved that the expresion in squared brackets is independent of s, say 
li{n). Applying the last equality to the orthonormal functions ■0n(s) and taking into 
account (NC3) and (NC4) we get 

^2n ^2n+2 



2n 2n + 2 
With the same technique we calculate 

L+(s, n - n) = /(s, n - l)g{s, n) + a{s) {/(s, n - 1) - g{s, n)} ^ 



Prom the properties between f{s,n) and g{s,n), the second term in the right side 
becomes zero. Substituting the values of these functions and H{s, n) we finally obtain 



L'^{s,n — 1)L {s,n) — 
+ n 



n 



'n-l 



n-l 



a{s) 



— a{s)H{s, n) 



It can be proved that the expresion in squared bracket is independent of s, say ^{n). 
Applying the last equality to the orthonormal functions ipnis) and taking into account 
(NC3) and (NC4) we get 



v[n] 



-Oln-\ln 



2n - 2 2n 

Obviously, v{n + 1) = /^(^t-)- These constants are given exphcitely in Table I 
Finally we have the desired relation equivalent to the Infeld-HuU-Miller factorization 
method: 

(NC5) L~{s, n + l)L+(s, n) = ^l{n) - cr{s)H{s, n) 

(NC6) L+(s, n)L-{s, n + 1) = /x(n) - a{s)H{s, n + 1) 

If we want L~^{s,n) and L~{s,n) mutually adjoint we have to divide both sides of 
(NC5) and (NC6) by 

^2n+2 ^2n 

2n + 2 2n 
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7. Factorization of difference equation of hypergeometric type 

For the case of orthonormal hypergeometric functions of discrete variable, we define 
from (NDl) the operator 



) = ^/{a{x) + T{x))a{x + l)E+ + ^/{a{x - 1) + t{x - 1)) 
- {2a{x) + t{x)) + Xn 
where E~^f{x) — f{x + 1), E~f{x) — f{x — 1), and the orthonormal functions satisfy 
H{x,n)(f)n{x) = 

As before we write the raising and lowering operators in the following way 

L~^{x,n) = u{x,n) + \/ {a{x — 1) + t{x — 1)) a{x)E~ 
L~{x, n) = v{x, n) + {(t{x) + t{x)) a(x + 1)E'^ 



where 



u{x, n) = a{x) 



XnTnix) ^ _^ _ ^^-j _ ^^^^ _ ^j-^-j 



n t'^ 2n 
Both expressions satisfy 

u{x + l,n) = v{x,n + 1) or 

u{x + 1, n — 1) = v{x, n) 

that can be proved by Taylor expansion. 
Now we calculate 

L~(x, n + l)L'^{x, n) — v{x, n + l)u{x, n) + {(t{x) + t{x)) a{x + 1) + u{x + 1, n) 
X \^^y{a(x) + T(x))a{x + l)E+ + ^y(a{x - 1) + t{x - 1)) a{x)E-^ 
Substituting the values for u{x,n),v{x,n) and H{x,n) we get 

L~{x, n + l)L"^(x, n) — 



~ ~f ~ ~/ '^V^ + ^) 



+ — ('^W + 



+ u{x + 1, n)H{x, n) 



It can be proved that the expresion in squared bracket is independent of x, say 
/i(n). Applying the last equality to the orthonormal function 0„(a;) and taking into 
account (NDl), (ND3) and (ND4) we get 

^^^^ = 2;r2;rr^"-^"+^ 
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With the same technique we calculate 

L^x, n - l)L'{x, n) = If.hl^Ini^ + _ 1 _ + A, 

X f-^^ + ^{x - Pn) + aix)) - {a{x - 1) + r{x - 1)) 



+ u{x, n — l)H{x, n) 



As before the expression in squared brackets is independent of x, say z/(n). Applying 
both sides of the last equality to the functions 0„(a;), and taking into account (NDl) 
(ND3) and (ND4) we obtain 

/ N ^2n-2 ^2n 

"^^^ = 2;rr2 2^7^"-^^" 
Obviously ^{n + 1) = /J-in) 

These constants are given explicitely in Table II. 

Finally the desired relations corresponding to the Spiridonov-Vinet-Zhedanov 
factorization method are 

(ND5) L"(x, n + l)L+{x, n) = ii{n) + u{x + 1, n)H{x, n) 

(ND6) L+{x, n)L~{x, n + 1) = fi{n) + u{x, n - l)H{x, n + 1) 

Again, if we want n) and L~{x,n) to be mutually adjoint, we have to divide 

both expressions (ND5) and (ND6)by 

^2n ^2n+2 

2n 2n + 2 

only in the case of Hahn and Chebyshev functions. 



8. Some comments 



The classical orthogonal polynomials we have presented in the preceeding sections are 
solutions of the second order differential equation 

(T{s)y';{s) + T{s)y'„{s) + A„?/„(s) = 

in the continuous case, or second order difference equation 

a{x)AVyn{x) + T(a;)Ay„(x) + A„y„(a;) = 

in the discrete case for uniform lattices where a{x) and t{x) are polynomials of at most 
the second and first degree respectively . 

Atakishiev and collaborators have generalized the classical orthogonal polynomials 
using a characterization based on the difference equation of hypergeometric type that 
covers all the cases defined by Andrews and Askey [16]. This characterization covers 
the q-analogue of classical orthogonal polynomials on non-uniform lattices. 
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Our paper should be implemented with the construction of raising and lowering 
operators for the orthogonal polynomials on non-uniform lattices in particular the q- 
analogue of the classical orthogonal polynomials. For this purpose we have at our 
disposal the analogue of difference equations, Rodrigues formula, recurrence relations 
for the orthogonal polynomials on non-uniform lattice, given explicetely by Nikiforov, 
Suslov and Uvarov [4]. 

An other approach for the same problem is given by Smirnov, via the factorization 
method suggested by Schrodinger for the solution of second order differential equation of 
hypergeometric type. Smirnov has applied this method to the finite difference equation 
on uniform lattices [9] and on non-uniform lattices [17], [18]. In his approach the raising 
and lowering operators are defined with respect to two indices: the first one, the degree 
of the orthogonal polynomials, the second one the order of the finite derivative with 
respect to the discrete variable. For this reason his raising and lowering operators are 
not equal to ours. 

A final comment to Tables I and II. In an unpublished Report of R. Koekoek 
and R.F. Swarttouw [20] tables are presented for orthogonal polynomials of the Askey- 
scheme and its q-analogue; among them one finds the raising and lowering operators 
of classical orthogonal polynomials of hypergeometric type. There are two points by 
which our tables are different from theirs. First, wc have calculated the raising and 
lowering operators from Rodrigues formula (sec (1, 2, 13) and (2, 2, 10) of Ref. [4]), 
but their raising and lowering operators are connected with some rccTirrcnce relations 
(see (1, 4, 5) and (2, 4, 13-17) of Ref. [4]) which are defined with respect to two indices. 
Besides that their tables do not cover the differential/difference equations, recurrence 
relations and raising/lowering operators with respect to the orthonormal functions of 
hyper geometric type as given in our tables. 

Acknowledgments 

The author wants to express his gratitude to Prof. Smirnov for valuable conversations 
and to Prof. A Ronveaux for his advise and encouragement and for bringing to his 
attention the papers of A. P. Magnus and G. Bangerezako. He is also very thankful 
to the Referees for their suggestions and new references. This work has been partially 
supported by D.G.I.C.Y.T. under contract #Pb96-0538 (Spain). 



14 



TABLE I 

Orthogonal polynomials of continuous VEiriable 



Hermite polynomials 
Hel) - 2sH'^{s) + 2nH^{s) = 

He2) sHn{s) = ^//n+i(s) + nH^^iis) 
He3) H^+,is) = 2sH^is)-H'^is) 
He4) H^_,{s) = ^H'^is) 

Laguerre polynomials 

Lai) sLfis) + (! + «- + nL^X-s) = 

La2) (n + l)L^+i(s) + (n + + (s - 2n - a - l)L^(s) = 

La3) (n + l)L^+i(s) = (s - n - a - 1)LJ:(5) + sL^ {s) 

La4) (n + a)L^i(5) = nL^{s) - sL^ {s) 

Legendre polynomials 

Lei) (1 - s')P;(s) - 2sP^(s) + n(n + l)P„(s) = 

Le3) (n + l)P„+i(s) = (n + l)sP„(s) - (1 - s')P^(s) 
Le4) nP„_i(s) = nsPn{s) + (1 - s2)P^(s) 

Jacobi polynomials 

Jl) (1 - 52)Pi"''^)"(s) + [/?-« - (a + + 2)5] Pi"'/^)'!.) 
+ n(n + Q; + /? + l)P^"''^)(s) = 

2(n+l)(n + a + /? + l) (.,^) 
^ (2n + Q; + /?+l)(2n + a + /3 + 2) ^ ^ 



+ 



+ 



2{n + a){n + (3) 



{2n + a + p){2n + a + p + 1) 

p^-a^ 

{2n + a + (3){2n + a + (3 + 2) 



P^"'/3)(s) = 



J3) 



2(ri + l)(7/ + n +/] + !) 



{2n + a + p + 2) 

' (n + a + P + l) 
_{2n + a + P + 2) 



{a- p) + {n + a + (3 + l)s 



J4) %^^^^^ef(.) 



(n + a + p + l) 



{2n + a + /3 + 2) 



{2n + a + p){2n + a + p + 2) 



{2n + a + /3) 

+ (2n + a + /3 + 1) - 

+ (l-s^)Pi"'^)(s) 
Normalized Hermite functions 

^n{s) = {2-n\^)-''\-'"hHr.{s) 

NHel) V'n(s) + (1 - s')V'n(s) + 2nV'n(s) = 

NHe2) V2(n + l)V'n+i(s) + A^V'n-i(s) - 2sV;„(s) = 

NHe3) L+(s, n)V'n(s) = #n(s) - = V2(n + 

NHe4) L-(s,n)^„(s) = s^/;„(s) + ^/'^(s) = v^t/'„_i(s) 

L~{s, n)L^{s, n)ipn{s) = 2{n + l)ipn{s) 



{P - a) - {n + a + P + l)s 



Normalized Laguerre functions 



T{n + a + l) 



NLal) silj'^is) + V^;(s) -^(^s + ^-2a-2^ V'„(s) + nV^„(s) = 

NLa2) y/{n + l){n + a + l)V'n+i(s) + + Q;)V'n-i(s) - (2n + a + 1 - s)V'n(s) 

1 



NLa3) L+(s,n)t/'„(,s) 



-(2n + a + 2 - s)Ms) - si^nis) 



= - V(^ + l)(n + Q; + l)V'n+i(s) 
NLa4) L'{s, nyipn{s) = --{2n + a - s)^/'n(s) + sij'^{s) = -^/n{n + a)'ipn-i{s) 



T{a + 1) 
1 



e-./2^a/2 



n-1 



V'o(s) = 

V^K«+l)n fe=o 

L'^(s, n - 1)-L~(s, n)'^n{s) = n(n + Q;)'0n(s) 
L~(s, n + n)'^n{s) = (n + l)(n + a + 



J]L+(s,n-l-A;)^o(s) 
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Normalized Legendre functions 



NLel) (1 - s'Xls) - 2s<(s) + n{n + l)V^n(s) = 

j2n + 1 /2n + 1 

NLe2) (n + 1) V ('''') + ^^^III^^-i^^) ~ + = 



/2r? + 1 

NLe3) L+{s,n)Ms) = {n + l)sMs) - (1 - = (n + 



/ 2n + 1 

NLe4) L-(s, n)V'n(s) = nsipn{s) + (1 - s^)V'n(s) = -^ ^n-l(g) 



^ n— 1 

V;„(s) = -^2^17 L+(s,n - 1 - k)Ms) 



fc=0 



L'^{s,n — 1)L {s,n)=n'^ 
L-{s, n + l)L+(s, n) = (n + 1)^ 



Normalized Jacobi functions 



^ M - , / ^'(^^ + + ^ + ^)(^ + + ^ + ^) (1 - g)f (1 + g)l (,,^) 

NJl) (1 - - 2.<(s) - i I - 2{a + + 

+ n{n + a + P + l)ip„{s) = 

2y/{n + l){n + a + l){n + p + l){n + a + p + l){2n + a + p + l) 
^ {2n + a + p + 2) V2n + a + /:/ + 3 ^"^^ ^^'^ 

2v/n(n + + + a + P){2n + a + P + 1) 
^ {2n + a + f3 + 2)^/2n + a + f3 - 1 'Vn-iis) + 

+ (2n + a + (3 + l)< ^^7- — — - s ) ^As) = 

^ ^ '\{2n + a + P){2n + a + P + 2) y""^' 

NJ3) L+(s, = {(n + a + /3 + 1). - - « - ri^) 

+ - a - (a + /3)s)|v'n(s) - (1 - = 

_ 2^{n + l)(n + g + l)(n + + l)(n + g + /3 + 1)(2to + g + /3 + 1) 

~ (2n + a + /3 + 2) V2n + a + + 3 V'n+ilsj 
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NJ4) L-{s, n)Us) = I - (n + a + /3 + 1)5 + - a - n') 
+ (2n + a + 3 + l) s - 

- i (/3 - a - (a + P)s) |Vn(s) + (1 - 

_ 2y/n{n + a){n + P){n + a + l3){2n + a + p + 1) 

~ {2n + a + /3) V2n + a + /3^ 'Vn-iis) 

% (s) = a + f3 + l ,^ _ ,| ,^ ,| 

° ^ ^ ^2"+/^+ir (a + 1) r (/3 + 1) ^ ^ ^ ^ 



I 2V(^ + l)(^ + a + l)(^ + /3 + l) 

^/{k + a + P + l)(2k + a + P + l) Jj^yJ 

r + / 1^r-^ ^/^^ 4.n(n + a){n + (3){n + a + P) 
L+{s,n-l)L {s,n)iJn{s)^ {2n + a + pf 



L {s,n + l)L^{s,n)ipn{s) 

_ 4(n + l)(n + g + l)(n + /3 + l)(n + g + /3 + 1) 
(2n + q; + /3 + 2)2 



TABLE II 

Orthogonal Polynomials of Discrete Vciriable 

Kravchuk polynomials 

Kl) ?!^L^kJx + 1) + XkJx - 1) + ^(P-g)-^^ feJa;) + I^fcja;) ^ 

n + 1 

K2) A;„+i(x) +p(A^ - n + l)A;„_i(x) + [n + p(A^ - 2n) - x] A;„(x) = 

Tl ~t~ 1 J? 

K3) kn+i{x) — -{x + n — N)kn{x) + xkn{x — 1) 

K4) p{N -n + l)^^{x + n- N)kn{x) + ^{N - x)kn{x + 1) 

Meixner polynomyals 

Ml) ijl{x + ^)mn{x + 1) + xmn{x - 1) - [nix + 7) + x] m„(x) + n{l — ij)mn{x) = 

M2) iimnj^i{x) — n{n + 7 — l)m„+i(a;) + [ii{x + n + 7) + n — a;] m„(a;) = 



M3) - //m„+i(x) = —piix + n + 7)m„(x) + xmn{x — 1) 

M4) — n{n + 7 — l)m„_i(a;) = —^{x + n + j)mn{x) + fj,{x + 'j)mn{x + 1) 

Charlier polynomials 

CI) IJ,Cn{x + 1) + XCn{x - 1) - {x + lj)Cn{x) + nCn{x) = 

C2) - ncn+i{x) - ncn-i{x) + (n + // - a;)c„(a;) = 

C3) - I^Cn+l{x) = -IXCn{x) + XC„(x - 1) 

C4) - ncn-i{x) = -jJiCnix) + Aic„(a; + 1) 



Chebyshev polynomials 

Tl) {x + l){N -X- l)tnix + 1) + a;(iV - x)tn{x - 1) - 

- [(AT - X - l)(x + 1) + x(7V - x)] tn{x) + n{n + l)tn{x) = 

T2) ^(n + l)tn+i{x) + {N^ - n^) tn-i{x) + ^(2n + 1) (iV - 1 - 2x) tnix) = 
T3) ^{n + l)tn+i{x) 



-{n + 1){N -2x -n-l)+ x{N - x) 



tnix) + x{N — x)tn{x 



T4) -n {N' - n') tn-i(x) 



hn + 1){N - 2x - n - 1) + n{n + 1) + {2n + 1) (x- ^ ^ 



- {x\\){N-x-\) 



tn{x) + {X+1){N-X- l)tn{x + 1) 



Hahn polynomials 

Hal) [x(N -X- (3-2) + ((3 + l)(N -1)] hl'^{x + 1) + x{N + a - x)hl'^{x - 1 
- [x{2N - 2x + a - /3 - 2) + (/3 + l)(Ar - 1)] /i^'^(x) 
+ n{n + a + P + l)hl'^{x)^Q 



Ha2) 



{n + l){n + a + l3 + l) p 
2n + a + l3 + 2 



^ (n + c.)(n + P)[N^n + a + /^)(iV - n)^.,, + (^^ + , + ^ + i) 



2n + a + /9 



X 



a - /3 + 27V - 2 



+ 



_ «2) (2Ar + a + /5) 
4(2n + a + /3)(2n + a + /3 + 2) 



x 



K'^{x) 
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Ha3) 



Ha4) 



(n + a + P + l) 



{P + 1){N -l)-{a + P + 2 + 2n)x 



2n + a + p + 2 
+ {N-n-(3-2)n + x{N + a - x)^h'^'^ 



{n + a){n + P){N + n + a + (3){N-n) p _ 



2n + a + /3 

x{N - X - /? - 2) + (/3 + 1)(7V - 1)1 + 1) 

+ 

+(A^ -n- (3- 2)n^ + n(n + a + + 1) + (2n + a + + 1) 

a-p + 2N -2 - a^) (2iV + a + /3) \ 



2n + a + ^\\ + 1)(/V - 1) - (. + ;3 + 2 + 2n). 



X 



4(2n + a + /3) (2n + a + /? + 2) . 



x(Ar - X - /3 - 2) + (/3 + 1)(A^ - 1) 



Normalized Kravchuk functions 



ln\{N-ny. / p^g 



kn{x) 



NKl) - ^)(^ + l)V'n(^ + 1) + ^^(N -X + l)x'ljjn(x - 1) + 

x{p-q)-Np n 

H ^ i^n{x) + -i^n{x) = 



NK2) jl{N-n){n+ l)iJn+i{x) + ^l{N-n+ l)n^„_i(a;) 
+^ [n + p(iV - 2n) - x] V'„(a;) = 



NK3) L+(x, = ^(x + n - N)il)n{x) + ^^(AT - x + l)xV'n(2; - 1) 



P 



(Ar-n)(n+l)V^„+i(x) 



NK4) L-(x, n)V'n(^) = ^(x + n - A^)V'n(a;) + J^{N - x)(x + l)V'n(a; + 1) 



= J^(iV-n + l)nV'„-i(x) 



x\{N-x)\ 



, n-1 



k=0 



P 

L'^{x, n — l)L~{x, n)ipn{x) = -{N — n + l)nijjn{x) 
L'~(x, n + l)L'^{x, n)'ijjn(x) = '^(N — n)(n + l)'0n(2;) 



Normalized Meixner functions 



NMl) ^l{x + -i){x + l)i)n{x + 1) + \/ ^IX{X + - l)^n{x - 1) - 

- [lJi{x + 7) + X] tpn{x) + n{l - IJ,)lpnix) = 

NM2) - VAi(n + 7)(n + l)iJn+i{x) - ^ ^n{n + 7 - l)^„_i(a;) + 
+ [//(x + n + 7) + n - x] -071(2^) = 



NM3) Ij^{x, n)ipn{x) = - [fi{x + n + 7)] ipn{x) + \/ lJix{x + 7 - IjV'nla^ ~ 1) = 
= - ^J n{n + -i){n + l)il)n+\{x) 

NM4) L-(x, n)V^„(x) = - [ii{x + n + 7)] ^n{x) + ji{x + l){x + '^)^ri{x + 1) 



= - •\/^(n + 7 - l)nV'n-i(a:) 



4)r 



L^{x, n — l)L~(x, n)'^n{x) = //(n + 7 - l)nV'n(2^) 
L~(x, n + l)L^{x,n)'^n{x) = iji{n + 'y){n+ l)ipn{x) 



Normalized Charlier functions 



Ipnix) 



ni V 



NCI) y/lJ,{x + l)lpnix + 1) + ^JWi^ni^ - 1) - + y^i'nW) + 'M'nWl = 

NC2) - \/ \x{n + l)V'„+i(x) - ^/]xh'4)n-\{x) + (n + - x)V'„(a;) = 



NC3) L+(a;,n)^/'„(a;) = -/iipnix) + ^/jlxi)n{x - 1) = ~\/ li{n + l)t/'„+i(a;) 



NC4) L (x, n)ipn{x) = -iJ.ip„{x) + s/ ii{x + + 1) = -^/Jmiln-l{x) 



ipo{x) 



XI 



i-iy 



n— 1 



:J]L+(x,n-l-A;)V'o(a^) 



L'^(x, n — l)L~(a;, n)V'n(2;) — iJ>nipn{x) 
L~{x, n + l)L"^(x, n)'^n{x) = //(n + l)'^n{x) 



Normalized Chebyshev functions 



NTl) {x + 1){N - X - 1) + yjnix + 1) + a;(iV - x)yjnix - 1) - 

- [(x + l)(Ar - x - 1) + x{N - x)] ipn{x) + n{n + l)ipn{x) = 



+ (2n + 1) ( ^^^^ - X ) = 



NT3) L+(x,n)ipn(x) ^ 



-{n + 1){N -2x-n-l)+ x{N - x) 



/.r N,/ ^ + 1 /(2n+ l)(iV2 -n2 -2n- 1) , 



NT4) L-(x,n)^„(x) = 



2n + 3 



-(n + 1)(A^ - 2x - n - 1) + n(n + 1) 
2 



TV — 1 

+ (2n + 1) ( X — ) - {x^\){N -x-\) 



^{X + l)(iV - X - 1)^^{X + 1) = f ""'^ ^n-M 
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n-l 



^ni^) = U{khy {2k + l){N^-l--2k-lf ^''^ n-l-k)\ M^) 



k=0 



71^ 



L^{x,n- l)L~{x,n)tpn{x) = —{N + n){N - n)ipn{x) 

in + 

L-(j;,n + l)L+(a;,n)V'n(a;) = ^ — -^{N + n + 1){N - n - l)'il)n{x) 



Normalized Hahn functions 



{2n + a + P + l)n\{N -n- l)\r{n + a + P + l) 
r{n + a + l)r(n + f3 + l)r{N + n + a + f3 + l) 



lr (N + a-x)r{x + P+l) 

r(7v-,T)r(,T + i) " ^' 

NHal) ^/{N -x- l){x + p + 1){N + a - x - l){x + l)^„{x + 1) + 

+ ^/{N - x){x + (3){N + a - x)xipri{x - I) - - x - l){x + P + 1) 
+ x{N + a-x) ^ipn{x) + n{n + a + l3 + l)ipnix) = 

y/{n + l){n + a + l){n + P + l){n + a + P + l)(2n + a + P + 1) 



NHa2) 



{2n + a + P + 2) 



{N + n + a + P + l)(N-n-l) 
""^1 2n + a + P + 3 ^"^^^"^^ 

^/n{n + a){n + P){n + a + P){2n + a + P + 1){N + n + a + P){N - n) 
{2n + a + P)y/2n + a + P - 1 
2N + a-p-2 



X ipn-i{x) + {2n + a + P + l) 



4 



iP^-a^)i2N + a + P) 

^ 4{2n + a + P){2n + a + P + 2) ^ ^^-^^'^ 



n + a + P+1 



NHa3) L+{x, n)iJn{x) = ^x{N + a - x){p + x){N - x)iJn{x - 1) - ^ ^ ^ ^ ^ 3 

X + 1)(A^ -l)-(a + P + 2 + 2n)x + (N - n - P - 2)n} 

,1 , , , J{n + l){n + a + l){n + P+l){n + a + P + l) 
+ x{N + a - x) iJn{x) ^ — /\ y /\ y I 



X 



(2n + a + /3 + 2) V2n + a + + 3 
V(2n + « + /? + l)(iV + n + a + /3 + l)(iV - n- l) 



NHa4) L-(x, n)V'n(2;) = V(x + 1)(A^ + a - x - + /3 + 1)(A^ - x - l)V'n(2: + 1) 
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+ 



n + a + (3 + l 



+ l)(iV - 1) - (a + /3 + 2 + 2n) X 



2n + a + (3 + 2 

+ {N -n- P- 2)n| + n{n + a + P + 1) + {2n + a + P + 1) 
2N + a-(3-2 {(3^ -a^){a + P + 2N) 



X 



4 4:{2n + a + p){2n + a + P + 2). 

_ (;V - . - 1)(. + /5 + 1)1 = ^!f±|f^^^ 

J {2n + a + p)\/2n + a + p — 1 

X V(2n + a + /3 + l){N + n + a + /3){N -n) ■ i{Jn-i{x) 



{a + (3 + l){N- l)\T{a + /?+!) T{N + a- x)T{x + /?+!) 



r{a + i)r{p + i)r{N + a + p + i) 



r{N-x)r{x + i) 



n-l 



k=0 



{2k + a + P + 2) 



^/{k + l){k + a + l){k + P + l){k + a + P + l) 



{2k + a + (3 + 3) 



{2k + a + P+l){N + k + a + (3+l){N -k-1) 
X L'^{x,n - 1 - A;) I • ■0o 

L'^{x, n — l)L~{x, n)il)n{x) 

_ n(n + a){n + (3){n + a + (3)(N + n + a + (3){N - n) 
~ (2n + a + (3y 

{n+l){n + a + l){n + (3 + 1) 



L {x,n + l)L'^{x,n)i(^n{x) 



{2n + a + (3 + 2)2 
X {n + a + (3 + 1){N + n + a + P + 1){N - n - l)ipn{x) 



References 

[1] Infeld L. and Hull T.E. "The Factorization Method" Rev. Mod. Phys. 23, 21-68 (1951). 

[2] Miller W. "Lie Theory and Difference Equations. 1" J. Math. Anal, and Appl. 28, 383-399 (1969); 

Lie Theory and special Functions, Academic, N.Y. 1968. 
[3] Nikiforov A.F. and Uvarov V.B. Special Functions of Mathematical Physics, Birkhauser, Basel 

1988. 

[4] Nikiforov A.F., Suslov S.K. and Uvarov V.B. Classical Orthogonal Polynomials of a Discrete 

Variable, Springer, Berlin 1991. 
[5] Atakishiev N.M. and Suslov S.K. "Difference analogs of the harmonic oscillator" Theor. Math. 

Phys. 85, 1055-1062 (1991). 
[6] Atakishiyev N.M. and Wolf K.B. "Aproximation on a finite set of points through Kravchuk 

functions" Rev. Mex Fis. 40, 366-377 (1994). 



24 



Atakishiev N.A., Jafarov E.I., Nagiyev S.M. and Wolf K. "Meixner oscillators" Rev. Mex Fis. 44, 
235-244 (1998). 

Atakishiyev N.M., Vicent L.E. and Wolf K.B. "Continuous vs discrete fractional Fourier transform" 

J. Comp. and Appl. Math. 107, 73-95 (1999). 
Smirnov Y.F. "On factorization and algebraization of difference equations of hypergeometric type" 
in Alfaro M. and al. Eds., Proceedings of the Int. Workshop on Orthogonal Polynomials in Math. 
Physics (On honour of Prof. Andre Ronveaux) Universidad Carlos III de Madrid, 1996. 
Bangerezako G., Magnus A. P. "The factorization method for the classical polynomials" in 
Priezzhev V.B. and Spiridonov V.P. Int. Workshop on Self-similar Systems J.I.N.R. Dubna 
1999. 

Bangerezako G. "The factorization method for the Askey- Wilson polynomials" J. Comp. and Appl. 

Math. 107, 219-233 (1999). 
Bangerezako G. Discrete factorization techniques for orthogonal polynomials on lattices 

(Dissertation) Un. Cath. Louvain, Louvain-la-Neuve, 1999. 
Lorente, M. "Quantum Mechanics on discrete space and time" in Ferrero M. and Van der Merwe 
A. eds. New Developments on Fundamental Problems in Quantum Physics^ (Kluwer, Academic 
Dordrecht 1997) pp. 213-224. 
Lorente, M. "Creation and annihilation operators for orthogonal polynomials of continuous and 

discrete variables" Electr. Trans. Num. Anal. (E.T.N.A.) 9, 102-111 (1999). 
Lorente, M. and Kramer P. "Representations of the discrete inhomogeneous Lorentz group and 

Dirac wave equation on the lattice" J. Phys A: Math. Gen 32, 2481-2497 (1999). 
Atakishiev N.M., Rahman, M. and Suslov, S.K. "On classical orthogonal polynomials", Constr. 

Approx.ll, 181-226 (1995). 
Smirnov, Y.F. "Finite difference equations and factorization method", Proc. 5th Wigner 

Symposium, Vienna, Austria (P. Kasperkovitz, D. Grau ed.) World Scientific, Singapure 1998. 
Smirnov, Y.F. "Factorization method: new aspects". Rev. Mex. Fis., 45, suppl. 2, pp. 1-6 (1999). 
Spiridonov, V. et al. J. Phys. A: Math. Gen. 30, L15-L21 (1997), Lett. Math. Phys. 29, 63 (1993), 
J. Phys. A 27, L669 (1994). 
[20] Roelof Koekoek, Rene F. Swarttouw "The Askey-scheme of hypergeometric orthogonal polynomials 
and its q-analogue" . Delft University of Technology, Faculty of Information Technology and 
Systems, Department of Technical Mathematics and Informatics, Report no. 98-17, 1998 (on 
line version at: _http://aw. twi.tudelft.nl/~koekoek/askey/index. html) 



